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ABSTRACT 

Our knowledge of the blazar surface densities and luminosity functions, which are fundamental 
parameters, relies still on samples at relatively high flux limits. As a result, our understanding of 
this rare class of active galactic nuclei is mostly based on relatively bright and intrinsically luminous 
sources. We present the radio number counts, evolutionary properties, and luminosity functions of 
the faintest blazar sample with basically complete (~ 95%) identifications. Based on the Deep X-ray 
Radio Blazar Survey (DXRBS), it includes 129 flat-spectrum radio quasars (FSRQ) and 24 BL Lacs 
down to a 5 GHz flux and power ^ 50 mJy and ~ 10^^ W/Hz, respectively, an order of magnitude 
improvement as compared to previously published (radio-selected) blazar samples. DXRBS FSRQ 
are seen to evolve strongly, up to redshift ~ 1.5, above which high-power sources show a decline 
in their comoving space density. DXRBS BL Lacs, on the other hand, do not evolve. High-energy 
(HBL) and low-energy (LBL) peaked BL Lacs share the same lack of cosmological evolution, which 
is at variance with some previous results. The observed luminosity functions are in good agreement 
with the predictions of unified schemes, with FSRQ getting close to their expected minimum power. 
Despite the fact that the large majority of our blazars are FSRQ, BL Lacs are intrinsically ^ 50 times 
more numerous. Finally, the relative numbers of HBL and LBL in the radio and X-ray bands are 
different from those predicted by the so-called " blazar sequence" and support a scenario in which HBL 
represent a small minority (« 10%) of all BL Lacs. 

Subject headings: galaxies: active — galaxies: evolution — BL Lacertae objects: general — quasars: 
general — radio continuum: galaxies — X-rays: galaxies 



1. INTRODUCTION 

Blazars are one of the most extreme classes of ac- 
tive galactic nuclei (AGN), distinguished by high lumi- 
nosity, rapid variability, high polarization, radio core- 
dominance (and therefore flat ^0.5] radio spectra), 
and apparent super luminal speeds. Their broad-band 
emission extends from the radio up to the gamma-rays, 
and is dominated by non-thermal radiation (synchrotron 
and inverse-Compton), likely emitted by a relativistic jet 
pointed close to our line of sight. This so-called "rela- 
tivistic beaming" gives rise to many interesting effects, 
which explain mo st blazar features (see the Appendix of 
lUrrv fc Padovanil [T9951V The properties of misdirected 
blazars are consistent with those of radio galaxies. In- 
deed, unified schemes fe.g.. lUrrv fc Padovanil [l995l ) as- 
cribe the differences between "beamed" (i.e., with their 
jets forming a small angle w.r.t. the line of sight) ob- 
jects and the so-called "parent population" to orienta- 
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tion effects. Within the blazar class, which includes flat- 
spectrum radio quasars (FSRQ) and BL Lacertae ob- 
jects, these are thought to be the beamed counterparts 
of high- and low- luminosity radio galaxies, respectively. 
The main difference between the two blazar classes lies 
in their emission lines, which are strong and quasar-like 
for FSRQ and weak or in some cases outright absent in 
BL Lacs. 

As a consequence of their peculiar orientation with re- 
spect to our line of sight, blazars represent a rare class of 
sources, niaking up considerably less than 5% of all AGN 
(|Padovanilll997f) . Therefore, previously available blazar 
samples suffered from small number statistics and rela- 
tively high limiting fluxes (until recently ~ 1 Jy and a few 
X 10~^^ erg cni^'^ s^^ in the radio and X-ray band respec- 
tively). The small size of these samples (^ 30 — 50 ob- 
jects) implies also that the derivation of beaming param- 
eters based on luminosity function stu dies, and therefore 
the viability of unified schemes (e.g., iPadovani fc Urrvl 
II99OI: lUrrv et all Il991f) is considerably uncertain, espe- 
cially at low powers. Moreover, as our understanding 
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of the blazar phenomenon is mostly based on relatively 
bright and intrinsically luminous sources, we have only 
been sampling the tip of the iceberg of the blazar popula- 
tion. For example, the best (and only!) radio luminosity 
function (LF) of BL L acs is still the 1 J y one, which 
is about 15 years old (jSticke l et alJ ll99lD. The situa- 
tion is only shghtly better for FSRQ. I Wall et all (|2005f ) 
have recently studied the FSRQ LF and evolut i on us - 
ing the Parkes 0.25 Jy sainple, \y hile iRicci et"al] (120061 ) 
had to use the iKiihr et al.l ()1981f ) 1 Jy sample to study 
the epoch-dependency o f the FSRQ LF. (The fact that 
lUrrv fc Padovanil (|1995( ) used a 2 Jy sample for their re- 
view paper gives an idea of how slow progress in this 
field is.) This state of affairs, and the lack of a size- 
able sample which includes both FSRQ and BL Lacs, 
has so far also prevented a specific test of the suggested 
possible evolutionary link between FSRQ and BL Lacs 
(|Cavaliere fc D'Enall200^ : iBottcher fc Dermej[200^ . 

The large majority of all known blazars have been 
discovered either in radio or X-ray surveys (but see 
[CoUing e et al. 2005, for a recent optically-selected BL 
Lac sample identified from the Sloan Digital Sky Sur- 
vey). Previous work has shown that X-ray and radio 
selection methods yield BL Lacs with somewhat differ- 
ent properties, especially as regards the frequency at 
which most of the synchrotr on power is emitted, z^poak 
(jPadovani fc Giommilll995"bh . The spectral energy dis- 
tributions (SEDs) of most radio selected BL Lacs peak, 
in a — vfi, notation, in the IR/optical bands, and these 
sources are now referr ed to as LBL (low-energy peaked 
BL Lacs: Padovani & Giommi| [r995b). By contrast, X- 
ray surveys select mostly HBL (high-energy peaked BL 
Lacs), whose energy output peaks in the UV/X-ray 
bands. The question of which of the two BL Lac sub- 
classes is the most numerous one has been a topic of 
debate in the past few years. This is not simply a de- 
mographical issue but touches upon the details of jet 
physics a nd the so-called "blazar sequence" (see below 
as well as lPadovanill2007l, for a rev iew). In this respect, 
we note that iPerlman et all (Il998f ) and iPadovani et alJ 
(j2002, 200i) have recently demonstrated that, contrary 
to previous (lack of) evidence and the predictions of the 
blazar sequence, FSRQ with SEDs similar to those of 
HBL, or HFSRQ, do indeed exist, although these sources 
fail to reach z^pcak values as extreme as those of HBL. 
One way to answer this population question is through 
relatively deep number counts. However, BL Lac num- 
ber counts are still based on samples put together in the 
early nineties and, therefore, at relatively high fluxes. 
Furthermore, while the strong evolution of FSRQ down 
to the (still relatively high) fluxes sampled seems es- 
tablished, the issue of BL Lac evolution is still quite 
open. X-ray selected samples (mostly HBL) have gen- 
erally shown indications of (small) negative evolution, 
i.e., with sources bei ng less lumin ous and /or less numer- 
ous i n the past (see iRector et al . 2000; B eckmann et al.l 
l2003l and refe rences therein; thi s result has been recently 
challenged bv ICaccianiga et all (2002 )). Radio-selected 
samples (mostly LBL), on the o ther hand, have ex hib- 
ited very weak, if any, evolution (jStickel et al.lll99i[) . It 
is important to notice that no sample, so far, has studied 
the evolutionary properties and LF of FSRQ, HBL, and 
LBL, within the same survey. 

Deeper, sizable blazar samples, have started to be 



assembled in the past f ew years by us ijPerlman et al.l 
119981 : iLandt et al. l I2001L hereafter P aper I and II re- 



spectively; see also Pa dovani et al.' ("2003^) and oth- 
ers (RGB: iLaurent-M uehleiscn et al. (199|h REX 
'Caccianiga et al." (2002); HRX: iBeckmann et aU (HODS'): 
CLASS: iCaccianiga fc Marchl (I2004D: and Sedentary 
iGiommi et al.l l2005l' see also " iPadovanil (|2002f ) for a re- 
view). Some of these samples take advantage of the 
fact that blazars are relatively strong radio and X-ray 
sources and therefore use a double radio/X-ray selection 
method^ . 

In Paper I and II we presented the methods, most of 
the identifications, and some preliminary results of the 
Deep X-Ray Radio Blazar Survey (DXRBS), a large-area 
(< 1,900 deg^ depending on X-ray flux) survey which 
reaches relatively faint X-ray (^ 2 x 10~^ erg cm~^ s^^) 
and radio (~ 50 mJy, depending on declination) fluxes. 
In this paper we present the radio number counts, evo- 
lutionary properties, and luminosity functions of blazars 
in the DXRBS sample, which is at present almost com- 
pletely (> 94%) identified. 

There are various features that make DXRBS a unique 
sample with which to address various issues: 

1. DXRBS is currently the faintest and largest blazar 
sample with nearly complete identifications; it 
therefore allows us to reach fainter and intrinsi- 
cally weak blazars, thereby providing further tests 
of unified schemes; 

2. DXRBS includes both FSRQ and BL Lacs within 
the same sample; it is then possible to compare 
the properties of the two classes, including num- 
ber densities, luminosity functions, and evolution, 
independently of selection effects due to different 
sample criteria and definitions. This is vital also 
to test the idea of an evolutionary link between the 
two classes; 

3. DXRBS includes both LBL and HBL within the 
same sample; this has never been achieved before 
and makes it possible to compare the properties 
of the two classes independently of obvious effects 
due to the different selection bands and methods, 
HBL being normally X-ray selected and LBL being 
typically radio-selected; 

4. DXRBS reaches radio (and X-ray) fluxes faint 
enough to provide d efinite tests to the so-called 
" blazar sequence" (|Fossati et al.l I1998D . which 
posits an inverse dependence of Vpeak on intrinsic 
power. This scenario not only envisions the non- 
existence of FSRQ with SEDs similar to those of 
HBL, which are present in DXRBS, but predicts 
also a dominance of HBL in X-ray and relatively 
deep radio surveys, which, until now, could not be 
proven or disputed due to the lack of deep samples; 

The DXRBS survey and its selection criteria, the sky 
coverage of the X-ray catalogue we have used, WGA- 
CAT, and the sample selection and identifications are 
described in § 2. In § 3 we analyze the DXRBS complete- 
ness, while § 4 describes the DXRBS number counts, and 

^ We put DXRBS and some of these surveys in perspective in 



DXRBS Counts, Evolution, and Luminosity Functions 



3 



§ 5 discusses the sample evolutionary properties. We ob- 
tain the DXRBS luminosity functions in § 6, while § 7 
summarizes our conclusions. 

Throughout this paper spectral indices are written 
5*1/ cx i/^" and the values Hq = 70 k m s~^ Mpc~^, 
r^M = = 0.3, and flA = 0.7 have been used (jSpergel et al.l 
|2003|) . To compare some of our results with previous 
work at times we have also adopted an Hq = 50 km 
Mpc~^, = 0, and JIa = (empty Universe) cosmol- 
ogy. Preliminary results on some of the topics addressed 
in this paper were presented by .Padovani (2001, . .2002,) . 

2. THE DXRBS SAMPLE 

The selection technique and identification procedures 
used for DXRBS have been described in Paper I and II. 
We summarize here our final sample selection, discuss the 
WGCAT sky coverage, and present our sample definition 
and identifications. 

2.1. Candidate Selection 

DXRBS takes advantage of the fact that all blazars are 
relatively strong X-ray and radio emitters. Selecting X- 
ray and radio sources with flat radio spectra (one of the 
defining properties of the blazar class) is therefore a very 
efficient way of finding these rare sources. By adopting 
a spectral index cut ar < 0.7 DXRBS selects all FSRQ 
(defined by < 0.5) and basically all BL Lacs, and 
excludes the large majority of radio galaxies. 

DXRBS initially was the result of a cross-correlation 
of all serendipitous (i.e., excluding targets) X-ray sources 
in the publicl y available j^Qg^T database WGA CAT95 
(first revision: IWhite. Giommi. fc Angelini 19951 ). having 
quality flag > 5 (to avoid problematic detections) , with a 
number of publicly available radio catalogues (the chosen 
X-ray and radio catalogues were selected on the basis of 
their large area and low flux limit). North of the celes- 
tial equator, we used the 6 and 20 cm G reen Bank sur- 
vey catalogues GB6 and N ORTH20CM (jGregorv et all 
Il996t IWhite fc Beckeijll992D . while south of the equator 
we used the 6 cm Parkcs-MIT-NRAO catalogue PMN 
(jGriflnth fc Wright 1993). For objects south of the ce- 
lestial equator, where a survey at a frequency different 
from the one of the PMN (6 cm) was missing when w e 
started this project (the NVSS [Condon et all (|1998l )]. 
now available, reaches in any case only 5 = —40°), we 
conducted a snapshot survey with the Australia Tele- 
scope Compact Array (ATCA) at 3.6 and 6 cm. This 
not only gave us arcsecond radio positions for our south- 
ern sources (we use the NVSS for the northern ones) but 
also radio spectral indices unaffected by variability. Note 
that the primary selection has been done at 6 cm, as the 
20 cm catalogues are used to derive spectral indices. 

A second version of WGACAT was released in May 
2000. In the course of the work on this version. White, 
Giommi, and Angelini became aware of a problem with 
the coordinates for 345 sequences in WGACAT95, caused 
by an error in converting the header of the event files. 
This error caused an offset in the source declination up 
to 1 arcmin but typically less. The fraction of sources 
affected is very small, ~ 1.4% and ~ 0.4% for the PMN 
and GB6 correlation respectively. For objects satisfying 
our completeness criteria (see § 12. 3p , no source was mis- 
takenly included because of this error, while only one had 
to be added. The Icr WGACAT positional errors (see Pa- 



per I for more details) range from 13 arcsec for the inner 
10' of the PSPC field to 53 arcsec for the 50 - 60' ring. 

The X-ray/radio matching was done as follows. WGA- 
CAT was correlated with both GB6 and PMN catalogues 
with a radius of 1.5 arcmin, excluding however sources for 
which the ratio between X-ray/radio offset and positional 
error was larger than 2 for the inner region (< 30 arcmin) 
of the iZ 05^ T" Position Sensitive Proportional Counter 
(PSPC) (as for such relatively large correlation radii one 
expects a non-negligible number of spurious matches; see 
§ 13. 2p . In the northern hemisphere the resulting sample 
was then correlated with the NORTH20CM catalogue 
with a radius of 3 arcmin, as the positional uncertainties 
of the NORTH20CM catalogue are considerably worse 
than those of the GB6 catalogue (160 arcsec at the 90% 
level compared to 10 — 15 arcsec at the la level respec- 
tively), and the 6 — 20 cm spectral index calculated. In 
the south we derived the 6—20 cm spectral index from 
PMN and NVSS data, summing up the flux from all 
NVSS sources within 3 arcmin from the PMN position 
for S > —40°, while we used our own ATCA observation 
to derive the 3.6 — 6 cm spectral index for S < —40° (this 
slightly different criterion translates in an estimated loss 
of only three blazar candidates: see details in Paper II 
and also § 12. 3[) . Finally, we excluded from our candidate 
list the following sources: clearly resolved out by our 
ATCA observations, with |6| < 10°, within 5° from the 
Large and SmaU Magellanic Clouds (LMC, SMC) and M 
31, and within 6° from the Orion Nebula. 

2.2. WGACAT Sky Coverage 

The sensitivity of the ROSAT PSPC instrument, be- 
sides the obvious dependence on exposure time and back- 
ground intensity, is a strong function of the position in 
the fleld of view. Consequently, the area of the sky cov- 
ered at any given flux (usually known as the sky cov- 
erage) is a complex function of flux. Two basic factors 
are responsible for the off-axis radius dependence: 1. the 
decrease of the instrument effective area at large offset 
angles (vignetting effect); 2. the degradation of the Point 
Spread Function (PSF) with distance from the center. 
Both factors can be described analytically by the fol- 
lowing relationships which give the minimum detectable 
count rate cr min in a PSPC image: 



cr„i„ = 900 VWi) r<10' (1) 

cr„,in = 900 y^ib/t) exp[(r - 10)/15] r > 10' (2) 

where t is the exposure time in seconds, r the off-axis 
radius expressed in arc minutes, and b the value of the 
local background (counts/pixel/s). The dependence on t 
and b, in background limited exposures {t > 2000 s), is 
given by the factor ^ (b/t) while the reduced sensitivity 
at large off-set angles due to vignetting and the PSF 
degradation is described by the exponential term exp[(r— 
10)/15] and is negligible at low (r < 10') off-axis angles. 

The form of eqs. ([T]) and ([2|) and the value of its pa- 
rameters have been derived by studying the dependence 
of WGACAT source count rates on t, b and r. Eqs. |T]) 
and 12]) are somewhat conservative since some sources 
can still be detected just below the threshold, but it en- 
sures that the number of spurious sources and source 
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confusion are reduced, so that WGACAT can be used 
for statistical studies. 

The WGACAT sky coverage has then been computed 
simply inverting the sensitivity laws ([1]) and ([2|) for each 
field. The following areas of the PSPC field of view have 
been excluded from the computation: 

1. r < 1.8', to exclude the target of the PSPC obser- 
vations; 

2. 13' < r < 24', to avoid the PSPC window structure 
which absorbs most of the photons in this region 
due to the wobble motion applied to most expo- 
sures; 

3. r > 45', to take into account the strongly reduced 
PSPC sensitivity at larger off-axis angles; 

4. 10% of the circular region between r — 24' and 
r = 45', to subtract the area covered by the eight 
equally spaced ribs extending outward from the in- 
ner ring. 

Finally, to take into account the DXRBS selection cri- 
teria (see § 13. ip . the sky coverage has been calculated 
excluding the WGACAT fields satisfying the following 
conditions: 

1. S> 75° 

2. -10° < 6 < 10° 

3. circular regions around M 31, LMC, SMC (i? = 5°) 
and the Orion Nebula (i? = 6°) 

4. (small) regions not covered by the PMN and the 
GB6 surveys. 

Count rates have been converted to 0.3 — 2 keV 
fluxes assuming a power law spectral model absorbed 
by an amount of neutral hydrogen equal to the Galac- 
tic value a s determined by the 21 cm measurements of 
iDickev fc Lockman (1990). 

Fig. □ shows the DXRBS/WGACAT sky coverage for 
various assumptions of the power law energy index, rang- 
ing from ttx = —0.2 to — 3.8, with a step equal to 0.4. 
Note how the area covered at a given flux depends on the 
assumed energy slope. To avoid the uncertainty intro- 
duced by assuming a single value, as often done, this was 
estimated individually for every source from the ha rdness 
ratios, as described in lPadovani fc Giommil (jl996D '^. 

2.3. Sample definition and identifications 

We have defined a complete sample from the sources 
which meet the following five criteria: 

1. ae-20 < 0.7 for S > -40° and as e-e < 0.7 for 
S < -40°;-* 

^ Our results are basically unchanged even assuming Qx = 1 or 
Ox = 1.5 for all our sources. 

This slight difference in the wavelength range used to derive 
the radio spectral index has very little effect on our sample. In 
Paper II, in fact, we showed that the mean difference between 
Q!3.6— 6 and 06—20 for core-dominated sources was 0.13. Expand- 
ing on the discussion done there, we point out that there arc no 
quasars with S < —40° and 0.5 < 03.6-6 < 0.63 satisfying our 
completeness criteria, so basically no FSRQ is lost. Moreover, the 
fraction of BL Lacs having Or > 0.7 is very small (~ 5% in the 
IPadovani. Giommi. fc Fiorj I I1997I) AGN catalogue). 




(0.3 - 2 keV) [c.g.s.] 



Fig. 1.— The DXRBS/WGACAT sky coverage for various as- 
sumptions of the power law energy index, ranging from Ox = —0.2 
(bottom left curve) to Ox = 3.8 (top left curve), with a step equal 
to 0.4. 

2. |6| > 10°, avoiding also the Large and Small Mag- 
ellanic Clouds, the Orion Nebula, and M 31; 

3. /20cm > 150 mJy for < (5 < 75°; 

4. /6cm > 51 mJy for -87.5° < ,5 < -37° and -29° < 
S <0°, /6cm > 72 mJy for -37° <S < -29°; 

5. 1.8' < PSPC offset < 13' and 24' < PSPC offset < 
45'. 

Point 1 has been discussed in §[231 point 2 excludes the 
Galactic plane and the usual bright, extended sources, 
points 3 and 4 will be addressed in § [3l while point 5 
is related to the derivation of the WGACAT sky cover- 
age (see § 12. 2p . These criteria were chosen in order to 
ensure that a well defined, flux-limited sample could be 
obtained. 

The above defining criteria were met by 219 blazar 
candidates, of which 77 were previously known sources. 
Details on the identification of the optical counterpart 
and spectroscopic observations for most sources, includ- 
ing objects not belonging to the complete sample, are 
given in Papers I and II. Identifications for the remain- 
ing objects and the final list of DXRBS sources will be 
presented in a future publication. 

We follow here a classification scheme slightly changed 
from that a dopted in Paper II, using the work of 
iLandt et all (|200al200l . The blazar class includes BL 
Lacertae objects, historically characterized by an almost 
complete lack of emission lines, and the flat-spectrum ra- 
dio quasars (FSRQ), which by definition display broad, 
strong emission lines. The separation of BL Lacs from 
radio galaxies, on the one hand, and from radio quasars, 
on the other hand, is somewhat controversial. Both BL 
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TABLE 1 

Complete sample composition. 



Class 


Total 


Newly 


Previously 


No redshift 






identified 


known 




FSRQ 


129 


79 


50 





BL Lacs 


24 


15 


9 


7 


SSRQ^ 


33 


24 


9 





Radio Galaxies'* 


17 


9 


8 





Unidentified [a^ > 0.5] 


16 [6] 


16[6] 




16[6] 



^ The SSRQ and radio galaxies samples are obviously incomplete, as our 
radio spectral index cut excludes by definition the majority of these 



Lacs and radio galaxies can have no or onl y narrow emis- 
sion li nes and therefore, along the lines of iMarcha et al.l 
we have used the value of the Ca break, a stel- 
lar absorption feature in the optical spectrum defined by 
C = (/+ — /_)//+ (where /_ and /+ are the fluxes in 
the rest frame wavelength regions 3750 — 3950 A and 
4050 — 4250 A respectively), to differentiate between the 
two. We have ad o pted a separation value of C = 0.4, 
since iLandt et alJ (|2GG2) showed that above this value 
sources become increasingly lobe-dominated (see their 
Fig. 6). 

We have adopted a dividing value of full-width half 
maximum FWHM — 1000 km/s between "narrow" and 
"broad" emission lines. We also make the commonly ac- 
cepted distinction between steep spectrum radio quasars 
(SSRQ) (a,. > 0.5) and FSRQ {a, < 0.5). In order to 
separate BL Lacs (with broad emission lines) from radio 
q uasars we have used the physical classification scheme 
of lLandt et al.l (|2004f) which separates sources into weak- 
lined and strong-lined radio-loud AGN based on their lo- 
cation in the rest frame equivalent width plane of the 
narrow emission lines [Oil] A3727 and [OIII] A5007 (see 
their Fig. 4). If our spectrum did not cover the posi- 
tions of both these emission lines we have adopted for 
BL Lacs (i.e., weak- lined radio- loud AGN) a limit of 5 
A on the rest frame equivalent width of any detected 
emissioi i line, which is the same as that used for the 1 Jy 
sample (IStickel et al.lll99lD . 

Tab. [T] gives the breakdown of the sample. One hun- 
dred and fifty-three sources turned out to be blazars, that 
is 129 FSRQ and 24 BL Lacs respectively. 

3. COMPLETENESS OF DXRBS 

The completeness of our sample is a complex function 
of many factors, namely: 

1. the completeness of the catalogues which make up 
DXRBS; this is clearly beyond our control but 
needs to be taken into account to the best of our 
knowledge; 

2. our cross-correlation radii; one here has to find 
a reasonable compromise between too large of a 
value, which produces a high number of spurious 
sources, and too small of a value, which increases 
incompleteness; 

3. our double X-ray/radio selection, which could re- 
sult in some incompleteness when comparing our 
results with purely radio-selected samples; 

4. our radio spectral index cut; 



5. the still missing identifications; 

6. a subtle second order effect, to do with non- 
serendipitous sources being included, which could 
result in an excess of sources at large radio fluxes. 

We discuss these points in turn below. 

3.1. Input catalogues 
3.1.1. GB6 

The GB6 catalogue (jGregorv et al.l I1996D covers the 
northern sky at 0° < S < +75° with a flux limit at 6 cm 
which is declination dependent and « 25 mJy. Since all 
our northern sources have /20cm > 150 mJy (see below), 
our spectral index cut of a,- < 0.7 implies /ecm > 65 
mJy, well above the GB6 limit. Indeed, all our northern 
sources have /ecm > 87 mJy. 

3.1.2. NORTH20 

The NORTH20 catalogue (jWhite fc Beckeilll992l ) cov- 
ers the sky at —5° < 5 < -1-82° with a formal flux limit at 
20 cm of 100 mJy. However, as described in the original 
paper, the catalogue is known to be incomplete below 
150 mJy. We checked on this by using the NVSS. Out 
of the 32,530 sources in the NVSS with flux > 100 mJy, 
-5° < (5 < -H82°, and \b\ > 10°, only ~ 72% have a 
NORTH20 counterpart. For fluxes > 150 and 300 mJy 
this fraction goes up to ~ 90% and ~ 97%, respectively, 
which shows that there is still some incompleteness above 
150 mJy. To remedy this we cross-correlated the WGA- 
CAT/GB6 sample with the NVSS using a radius of 1.5 
arcmin. Spectral indices were derived as described in 
Padovani et al. (2003), that is by summing up the 1.4 
GHz flux from all NVSS sources within a 3 arcmin radius 
(corresponding roughly to the beam size of the GB6 sur- 
vey). This resulted in eleven sources with /20cm > 150 
mJy, PSPC offset outside the 13 and 24 arcmin range 
(ii l2.3p . and ar < 0.7 being added to our sample. We can 
then safely consider DXRBS complete for /20cm > 150 
mJy and 0° < S < +75° (the region of the sky covered 
by both GB6 and NORTH20 catalogues). 

3.1.3. PMN 

The PMN catalogue (jGriffith fc Wright! [l99l ) com- 
prises four different surveys in the —87.5° < S < +10° 
range, with flux limits at 6 cm which are in most cases 
declination dependent and extend down to ^ 20 mJy. As 
the GB6 catalogue reaches S = 0°,we did not include any 
PMN "northern" source. Also, for simplicity we included 
in our cross-correlation only sources with /gem > 51 mJy 
in the area covered by the Southern, Tropical, and Equa- 
torial surveys (-87.5° < 6 < -37° and -29° < 6 < 
+ 10°), which is the maa;zmw77i declination-dependent flux 
limit for these surveys, and sources in the Zenith survey 
(-37° < S < -29°) with /gem > 72 mJy, which is its 
limit. Apart from the declinations covered by the Zenith 
survey, we are then well above the formal completeness 
limits of the PMN survey. 

3.1.4. WGACAT 

The first versio n of WGACAT 

(jWhite. GiommiT M Angclini! [T995h covers ^ 10% 
of the sky to varying degrees of sensitivity, with ex- 
posures typically a factor of 100 longer than those 
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achieved during the six month ROSAT All Sky Survey 
(RASS). Since we cross-correlate WGACAT with radio 
catalogues, the chance that one of our sources is a 
spurious X-ray detection is vanishingly small. The 
DXRBS X-ray flux limits depend on the exposure time 
and the distance from the center of the PSPC (see 12. 2p 
and vary between ~ 10^^"' and ^ 10^^^ erg cni^^ s^^. 

3.2. Cross-correlations 

A cross-correlation between two catalogues will pro- 
duce a number of spurious associations Ngp which de- 
pends on the number of sources in the first (X-ray) cat- 
alogue Ni, the surface density in the second (radio) cat- 
alogue 5*, and the correlation radius r as follows 

Nsp = 8.73 X 10-^ (r/arcmin)^ {S/[#/deg'^]) Ni. (3) 

The number of spurious matches in a given thin shell 
of radius Ar will be 

A^sp/Ar 1.75xl0"3(r/arcmm) {S / [# / deg'^]) Ni. (4) 

The completeness of our cross-correlation was com- 
puted in a way which is as independent as possible of "a 
priori" assumptions, including the positional errors, but 
relies only on the previous equations. Namely, the cross- 
correlations were done with a radius larger than the one 
we finally adopted (typically 5 arcmin) and the matches 
were binned in radial shells (typically 6 arcsec) as N/Ar. 
The results were then plotted as a function of the posi- 
tional offset between the matches. This is shown for the 
WGACAT - GB6 correlation in Fig. [U (top), which dis- 
plays the strong peak at small radii due to real matches, 
followed by a decline, and then by the rise oc r due to 
the predominance of spurious sources at large radii. 

The total number of spurious sources for a given corre- 
lation radius could be estimated using eq. However, 
this requires a reliable estimate of S, the surface den- 
sity of sources in the radio catalogue, which might not 
be straightforward to determine when the flux limit de- 
pends on position, as is the case for both GB6 and PMN 
catalogues. Most importantly, the coefhcients in eq. ^ 
and (|4]) neglect clustering, which would slightly increase 
Nsp. Therefore, we left the numerical coefficient free and 
found the best-fit value by imposing that above a large 
enough radius (determined interactively) all sources are 
spurious, that is Ngp/Ar = N/Ar (solid line in Fig. [2] 
[top]). This gives us also an estimate of the complete- 
ness level at every correlation radius, which is defined as 
A^roai(r)/Afroai(r ^ oo), where Ar,eai(r-) = N{r) - Nsp{r) 
(in practice A^rcai('') converges at r > 3 arcmin). 

Fig. [2] (bottom) shows the dependence of completeness 
and real matches (1 — spurious) fraction as a function 
of positional offset for the WGACAT - GB6 correlation. 
We adopted as our correlation radius the value at which 
the two curves cross, that is 1.5 arcmin (the same result 
apphes to the WGACAT - PMN correlation). That is, 
we are maximizing our completeness while at the same 
time minimizing the fraction of spurious sources. Our 
results are as follows: the completeness level is ~ 87% 
for the WGACAT - GB6 correlation and ~ 85% for the 
WGACAT - PMN correlation. The percentages of spu- 
rious matches are ^ 13% and ~ 17% for the GB6 and 
PMN correlations respectively, very close to the fraction 
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Fig. 2. — Top: the number of matches per arcmin shell as a func- 
tion of positional offset for the WGACAT - GB6 correlation. The 
strong peak at small radii is due to real matches, while the linear 
increase at large radii is due to spurious sources. The solid line 
represents the number of spurious sources as predicted by eq. (|4]l, 
fitted by assuming that at large radii all sources are spurious. Bot- 
tom: the dependence of completeness (solid line) and real matches 
(1 — spurious) fraction (dashed line) on offset (see text for details). 
The dotted line at 1.5 arcmin defines our adopted cross-correlation 
radius. 

of sources with ratio between X-ray/radio offset and po- 
sitional error larger than 2 13% and ~ 21%), which 
were already excluded from the sample (§ 12. 1[) . For our 
adopted cross-correlation radius of 3 arcmin, the com- 
pleteness level of the WGACAT/GB6 ~ NORTH20CM 
correlation is ~ 99% (this very high value is explained 
by the fact that we are correlating two radio catalogues 
at similar frequencies). As discussed above, the posi- 
tional uncertainties of WGACAT sources depend on their 
PSPC offset, so one might worry that in the outer parts 
we could be more incomplete (see a preliminary discus- 
sion in Paper I). However, the completeness level for the 
WGACAT - PMN correlation for PSPC offsets between 
30 and 45 arcmin is ~ 83%, not significantly different 
from that of the fuU WGACAT (~ 85%). 

To take the incompleteness related to the adopted 
cross-correlation radius into account all our number 
counts and luminosity functions have been multiplied by 
the factor 1/0.85 - 1.18. 

3.3. The double X-ray/radio selection 

As many recent blazar surveys, DXRBS is not simply 
radio or X-ray flux-limited, like the previous "classical" 
blazar samples, but a dopts a do u ble radio/X -ray selec- 
tion. As discussed bv lPadovtmil ()200lL [20M) . it is im- 
portant to assess what region of parameter space such 
a survey is sensitive to, in order to understand what 
constraints it can or cannot put on blazar demograph- 
ics. In particular, one should ask if such a sample can 
provide meaningful radio or X-ray constraints (number 
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log (0.1 - 2.4 keV) [erg cm ^ s '] 

Fig. 3. — Top: The sampling of the radio flux - X-ray flux plane 
by different blazar surveys. Thick lines represent "hard" survey 
limits, while thin lines are the fluxes reached in a band other than 
the one of selection. Sources belonging to a given survey occupy 
a region of the plane whose bottom-left corner is indicated by the 
thick/thin lines. The long-dashed line at /x//r = 3.2 X 10~ erg 
cm~^ s~^ Jy~^ separates sources with cti-x < 0.78 (high-energy 
peaked blazars) from sources with ctrx > 0.78 (low-energy peaked 
blazars), while the hatched regions represent the "forbidden" zones, 
different for BL Lacs and FSRQ at the high /x//r end, where no 
blazars have been found so far. Bottom: The short-dashed lines 
show the parameter space sampled by DXRBS. The diagonal solid 
line defines a DXRBS subsample with complete coverage of the 
plane for fx/ fv > 4 X 10~^^ erg cm~^ s~^ Jy"'^, while the vertical 
dot-dashed line defines a complete X-ray fiux limited DXRBS LBL 
sample with /x > 1.6 X 10~ erg cm"'^ s~^. See text for more 
details. 

counts, luminosity functions, etc.), i.e., if it can provide 
a representative blazar sample which can be considered 
flux-limited in one band. The answer to this lies in how 
a survey covers the blazar radio flux - X-ray flux plane 
and how this compares to the region of this plane which 
can be occupied by blazars. 

To this aim, we have estimated the range of /x//r 
covered by blazars using the multiwavelength AGN 
catalo gue put together by iPadovani. Giommi. &: Fiori 
(Il997l). add i ng ad ditional X-ray information from 
iSiebert et al.l ()1998[ ). We get a value for the smallest 
X-ray-to-radio flux ratio 8 x 10~^^ erg cm~^ s~^ Jy~^j 
for both FSRQ and BL Lacs (upper diagonal dotted line 
in Fig. [3]). As regards the largest value, this differs 
by about one order of magnitud e for the two classes, a s 
expected based on the results of IPadovani et al.| ()2003D . 
being fx/ fr ~ 2 x 10~^ erg cm~^ s~^ Jy~^ for BL Lacs 
and /x//r ~ 3 X 10-1° erg cm'^ Jy^^ for FSRQ 
(lower diagonal dott ed lines in Fig . [H). T his agrees with 
the values found bv lGiommi et al.1 (|2006f) by correlating 
the NVSS sources with > 1 Jy with the RASS. 

Thes e i? 05*^ T* base d results hold even at fainter X-ray 
fluxes. iBassett et al.l (|2004.1 present X-ray data for 15 
high-redshift radio-loud quasars, 12 of which have radio 



spectral index information and classify as FSRQ (this is 
by no means a complete sample, although it represents 
~ 30% of all currently known z > 4 radio-loud quasars). 
The minimum /x//r value for this sample is ^ 10"^^ 
erg cm"^ s"^ Jy~^, i-e., basically the same as the one 
obtained based on ROSAT data, despite the fact that 
the typical X.-T: ay flux for this sample is ~ 10"^^ erg cm~^ 
s~i, i.e., about one order of magnitude smaller than for 
DXRBS. 

Fig. |3F (top) shows the sampling of the radio flux 
- X-ray flux plane by DXRBS compared to that b y the 
"classical" blazar surveys (1 Jy (jStickel et al1[l99ll).\Em- 
stein Medium Sens i tivity Survey (EMSS) (|Stocke et all 
119911: iRector et al.l l200Clt) , and E instein Imagiiig Pro - 
portional Counter (IPC) Slew (jPerlman et all Il996f) 
BL Lac samples and the 2 Jy FSRQ samp le 
(iWaU fc Peacockl [l985l : Idi Sereeo Ahghieri et al.l 11994)) 
and a number of recent surveys with double radio/X- 
ray selection (RASS - Green Bank (RGB) survey 
(jLaurent-Muehl eiscn et al. 19991 , the Radio-Emitting X- 
ray (REX) survey (Caccianiga ct aL 199ijj 2002), and the 
Sedentary survev (iGiommi. Meima. fc Padovanil 119991 : 
iGiommi et al.ll2005l )). Note that all of these recent sur- 
veys are primarily devoted to BL Lacs, although REX 
inclu des also emission line sources and IPadovani et al.l 
(1200^) have cross-correlated the RGB sample with the 
NVSS to obtain radio spectral indices and extract the 
FSRQ. We also add a very recent FSRQ survey which, al- 
though not as deep as DXRBS and purely radio-selected, 
goes deeper than the 2 Jy sample, namely the Parkes 0.25 
Jy flat-spectrum sample (Wall ct al. 2005). To the best 
of our knowledge, all other recent surveys which reach 
fainter radio fluxes than DXRBS are either not com- 
pletely identified and/or have a relatively bright optical 
magnitude cut which translates into severe constraints 
on the number and type of sources included. 

In Fig. [3] (top) every survey is characterized by ei- 
ther one or two flux limits (thick lines), while the small- 
est flux reached by a sample in a band other than the 
one of selection is given by a thin line. For example, 
the EMSS reaches /x ~ 2 x IQ-^^ erg cm'^ s'^ (thick 
line), by default the X-ray selection limit (actually, the 
faintest of various X-ray limits, due to the serendipitous 
nature of the survey). A limit in one band translates 
into a limit in the other one and in this case the ra- 
dio faintest EMSS BL Lac has a flux /r ~ 1 mJy (thin 
line). The sources of a given survey occupy a region of 
the flux-flux plane whose bottom-left corner is shown in 
the figure. The long-dashed line in the figure (X-ray- 
to-radio flux ratio /x//r = lO"^^'^ erg cm~^ s"^ Jy~^ 
or Oi-x ^ 0.78) divides high-energy peaked from low- 
energy pe aked blaza rs (bo th BL Lacs and FSRQ; see 
IPadovani et al. 200 [20031 ) 

Fig. |3| clearly illustrates the limitations of surveys with 
double radio/X-ray flux limits. A survey whose limits 
fall quite far from both dotted lines will not provide a 
complete picture of the blazar population. For example, 
the REX survey cannot provide BL Lac radio number 
counts to be compared with the predictions of a beaming 
model based on the 1 Jy sample, simply because it does 
not include all the BL Lacs above its radio limit (as it 

^ This is an updated and re vised version of a figure originally 
shown in lPadovanH 1I2OOII. l200l) . 
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misses all those above the radio limit but below the X-ray 
limit). For the complementary reason, neither can REX 
provide X-ray number counts to be compared with the 
predictions from a beaming model based on the EMSS 
sample. REX will provide radio number counts for HBL, 
given its proximity to the HBL/LBL dividing line, and 
X-ray number counts for LBL (as it detects all LBL above 
its X-ray flux limit). The same arguments apply to the 
RGB survey, which has the further problem of an optical 
limit {B < 18). In short, surveys with two flux limits can 
provide a complete picture of the blazar population only 
if one of the limits is relatively close to the edge of the 
region of parameter space occupied by blazars. 

DXRBS misses sources only in a small corner of the 
radio flux - X-ray flux plane between its X-ray flux 
limit and the lower limiting X-ray-to-radio flux ratio for 
blazars (Fig. [31 bottom). For a value of this parameter 
'--^ 8 X lO"^"* erg cm~^ s"^ Jy~^ and a radio limit of 51 
mJy, in fact, DXRBS should have been sensitive to X-ray 
fluxes down to ~ 4 x 10~^^ erg cm~^ s~^ to be consid- 
ered equivalent to a purely radio flux-limited sample. It 
then follows that no source with 51 < /ecm < 250 mJy 
and /x < 2 X 10"^'' erg cm~^ s~^ (see Fig. |3]) could have 
been included in the survey. In fact, DXRBS covers fully 
the blazar region of the radio flux - X-ray flux plane 
at radio fluxes > 250 mJy, as for this radio limit and 
given DXRBS' X-ray limit, our survey touches the "for- 
bidden" zone in Fig. [3] (and happens to overlap with the 
PKS 0.25 Jy). Furthermore, DXRBS includes blazars 
with /r > 51 mJy and /x//r > 4 x 10"^'^ erg cm~^ s~^ 
Jy~^, as shown by the solid line in Fig. [31 bottom. We 
will address this point again in the following sections, but 
we stress here that the missed small corner of the radio 
flux - X-ray flux plane has no influence on the number 
counts, evolution, and luminosity function of DXRBS, 
since to evaluate these we properly take into account our 
flux limits. However, a correction is required when com- 
paring these properties to those of purely radio-selected 
samples (see § [4]). Finally, an X-ray flux limited LBL 
sample with /x > 1.6 x 10~^^ erg cm~^ s~^ (dot-dashed 
line in Fig. [31 bottom) can be defined within DXRBS 
(see Section |42)) . 

3.4. Radio spectral index cut 

Our sources are selected to have < 0.7 (ii l2.3p . How- 
ever, the radio data used to derive the ra dio spectral 
index are simultaneous only for S < —40°. IWall et aTl 
()2005() have recently addressed this issue and pointed 
out that any flux-limited survey will preferentially se- 
lect sources in an up-state, whereas flux density mea- 
surements taken at a different time reflect sources in a 
mean state. They also point out that it is the variations 
at frequencies above the survey frequency which matter. 
Since our radio spectral indices have been derived using 
frequencies below the survey frequency, this would result 
in an artificial flattening of the spectra, with the inclu- 
sion of some extra sources in our sample. However, the 
fact that all our data come from flux-limited surveys is 
bound to mitigate this effect. 

3.5. Missing identifications 

Our sample is not yet fully identified. At the time of 
writing (January 2007), there are still sixteen sources in 



the complete sample for which we do not have a classi- 
fication (seven in the southern and nine in the northern 
hemisphere). Out of the 16 still unidentified objects, 
we note that six have ar > 0.5 and therefore cannot 
be FSRQ. Based on the relative fraction of the different 
DXRBS classes, we estimate that we are missing « 9 
FSRQ and w 2 BL Lacs, with the remaining five sources 
being SSRQ (4) and radio-galaxies (1). Therefore, we 
are likely to stiU be missing w 7% and « 8% of DXRBS 
FSRQ and BL Lacs respectively. We address this incom- 
pleteness in § [5] 

3.6. Non-serendipitous sources 

Although we have obviously excluded all i?OS'^T tar- 
gets from our samples, a subtle, second-order effect is 
present which results in the inclusion of additional non- 
serendipitous sources. Consider the case of an X-ray tar- 
get which was discovered because in the field of view of 
a well-known and, therefore, likely radio-bright blazar. 
When this source is observed as a ROSAT target, the 
quasar will appear to be serendipitously in its field of 
view, while in fact it is not there by chance and needs to 
be excluded. We have then done so for the easiest cases 
(e.g., radio- loud objects in fields of Einstein Medium Sen- 
sitivity Survey sources), also by going back to the original 
i? 05^ T" proposal, but in some other cases it is very dif- 
ficult to ascertain the relation between the target and 
a radio source in the field. This might translate into a 
residual excess at high radio fluxes (see § [4]) . 

4. NUMBER COUNTS 

The area of the sky over which every source could be 
found, used for the number counts, volume calculations, 
and LFs, is determined by its X-ray and radio fluxes as 
follows: 

1- if /ecm > 72 mJy and /20cm > 150 mJy, the radio 
source is detectable over the whole sky; we then use 
the full WGACAT sky coverage and determine the 
appropriate area based on the values of the X-ray 
flux and X-ray spectral index; 

2. if 51 < /ecm < 72 mJy and /20cm > 150 niJy, the 
radio source is detectable over the whole sky ex- 
cluding the PMN Zenith area (-37° < <5 < -29°); 
we then use the WGACAT sky coverage for the 
whole sky excluding that region; 

3. if /ecm > 72 mJy and /20cm < 150 mJy, the ra- 
dio source is not detectable in the northern hemi- 
sphere; we then use the WGACAT sky coverage for 
the southern sky; 

4. if 51 < /ecm < 72 mJy and /20cm < 150 mJy, the 
radio source is detectable over the whole southern 
sky excluding the PMN Zenith area (-37° < S < 
-29°); we then use the WGACAT sky coverage for 
the southern sky excluding that region. 

As discussed in § 13. 3| DXRBS is " almost" equivalent 
to a radio flux-limited sample. In fact, it misses some 
(but not all) sources with 51 < /ecm < 250 mJy and 
/x//r < 4 X 10~^^ erg cm~2 s"^ Jy"^ (see Fig. [3|). 
We corrected the number counts for this effect by first 
deriving the /x/ /r distribution for the sub-sample with 
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Fig. 4. — The integral number counts for DXRBS BL Lacs at 
5 GHz (6 cm). The solid line represents the total number counts, 
while filled circles show the values at selected fluxes to show the 
errors involved. The dotted-dashed line shows the counts without 
the correction which we ap plied to make DXRBS equivalent to a 
radio flux-limited sample 13.31 1. The open square represents the 
surface density for the 1 Jy sample, while the open triangle repre- 
sents the surface density for the S5 sample. Both of these points 
have been corrected to take into account our somewhat different 
deflnition of a BL Lac. The dashed line and filled triangles are the 
number counts for the BL Lacs with «rx < 0.78, while the dotted 
line and flUed squares are the number counts for the BL Lacs in 
the HBL box. See text for details. 

/ecm > 250 mJy (weigh ted for the effect of the sky cov- 
erage; see Pa dovani et al . 2003). This sub-sample cov- 
ers the full range of /x//r for blazars and therefore its 
/x//r distribution should be unbiased. We then com- 
pared that to the fx/ fr distribution for the whole sam- 
ple and estimated the fraction of missed sources with 
/x//r < 4 X 10~^^ erg cm~^ s~^ Jy~^ as a function of 
/x//r- The /x//r values were then converted to radio 
fluxes by using our X-ray limit and the number counts 
were then corrected. (For example, the correction for 
/ecm > 100 mJy was derived by evaluating the fraction 
of additional sources with /x//r < 2x 10^^'^ erg cm~^ s^^ 
Jy~^, this being equal to our X-ray flux limit of 2 x 10~^^ 
erg cm~^ s~^ divided by this radio flux.) The correction 
is not large, being ~ 13% for 100 < /gem < 250 mJy 
and ~ 70% only for /ecm < 100 mJy, and within « la 
from the uncorrected values (compare the dotted-dashed 
and solid lines in Fig. [3]and[S]). Given the small number 
statistics for the BL Lac sample (only 8 BL Lacs have 
/ecm > 250 mJy), we evaluated this for FSRQ and then 
applied the same correction to the BL Lac sample. 

4.1. BL Lacs 

Fig. Ulpresents the integral number counts at 5 GHz for 
the DXRBS BL Lacs (s olid line) , compared to the values 
derived from the 1 Jy (iStickel et al.l ll99lL open square) 
and S 5 (jKiihr et al.l Il987l as updated by IStickel et al.l 
(|1993f) : open triangle) samples. The filled circles show 



the values at selected bins to show the errors involved 
without crowding the plot. In making the comparison 
with previous BL Lac surveys we note that their defini- 
tion of BL Lac was somewhat different from ours, which 
is less restrictive and reaches a,- = 0.7 (while previous ra- 
dio samples were defined by < 0.5). Since 2/3 (16/24) 
of the DXRBS BL Lacs also fulfill the 1 Jy definition, the 
1 Jy and S5 points were multiplied by the factor 3/2 to 
take this difference into account. 

A couple of points need to be considered before dis- 
cussing our results: 1. we do not characterize the un- 
certainties by simply considering the total number of 
sources, as the error bars would become progressively 
smaller, which is misleading. At relatively faint X-ray 
(and, on average, radio) fluxes, in fact, the surveyed area 
decreases and the number of sources gets smaller. We 
take that into account by summing in quadrature the 
errors for the individual sources. As a result, the error 
bars are larger at high and low fluxes, where the number 
of sources is small, and smaller at intermediate fluxes, 
where the statistics is better; 2. the radio flux limit is 
higher at 20 cm than at 6 cm. This has an impact on the 
northern counts only below 150 mJy, our adopted com- 
pleteness limit at 20 cm, as for ai- ^ (the median for 
our sample), /ecm ~ /20cm- Given the lower flux limit of 
the 6 cm (PMN) survey, the southern counts go deeper 
than the northern ones. 

Taking all this into account, we can say that our counts 
agree with previous estimates at relatively high fluxes 
(250 mJy). The BL Lac with the largest radio flux in 
our sample has /gem ^ 0.99 Jy so we cannot really com- 
pare our results with the 1 Jy sample. Given the fact 
that for the northern and southern hemisphere we have 
used different radio catalogues and a somewhat different 
selection, we have also checked that the southern and 
northern counts agree within the errors down to ~ 150 
mJy; below this value only southern sources contribute. 

4.2. The relative numbers of HBL and LBL 

Until about ten years ago HBL (then called XBL for 
X-ray selected BL Lacs) were regarded to have their jets 
seen at larger angles w.r.t. to the line of sight as com- 
pared to LBL ( then called RBL for ra dio-selected BL 
Lacs; see, e.g.. lUrrv &: Padovanil I1995L and references 
therein). This meant also that XBL/HBL were thought 
to be the most num erous subclass. In the mid 1990s 
'Gio mmi fc Padovanil (19941 ) and iPadovani fc Giommil 
y,995a') proposed a new interpretation for the existence 
of HBL/XBL. Against the prevailing view at the time, 
HBL/XBL were instead suggested to be intrinsically rare 
and to represent the 10% of BL Lacs with relatively 
high i^pcak- The fact that they were the dominant class 
in the X-ray band was thought to be a simple selection 
effect related to their SED. Namely, X-ray surveys were 
sampling the BL Lac radio population at relatively low 
radio fluxes and mostly detected the small fraction of 
objects with high /x/ /r ratio s. This scena r io is n ow also 
supported by the results of iLandt et al.l ()2002f ). which 
show no difference in jet orientation betwee n the two BL 
Lac subclasses (see also lRector et al.ll2000D . 

By the end of the 1990s the "different orientation" 
paradigm was replaced by a new scenario. The so-called 
"blazar sequence" proposed that LBL and HBL sampled 
the higher and lower part of the jet bolometric luminosity 
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function, respectively, i nstead of different p arts of the ra- 
dio luminosity function ([^ss^reF^iIl993). The "blazar 
sequence" , wh ich was later expan ded to incorporate also 
radio quasars (jFossati et al.|[l998l ). advocates once more 
that HBL are more numerous than LBL. Furthermore, 
since it is based on the assumption that an inverse de- 
pendence exists between Vpcak and intrinsic power due 
to the effects of the more severe electron cooling in more 
powerful sources, it excludes the existence of blazars with 
high radio powers and high Vpcak (i-e., HFSRQ). 

The predictions for the relative numbers of HBL and 
LBL in radio and X-ray surveys are dramatically differ- 
ent in the two scenarios and can be put to test with deep 
blazar surveys. In the first case (jPadovani fc Giommil 
119953*). HBL represent a minority, and their fraction in 
radio surveys should be constant and w 10%. More- 
over, X-ray surveys should detect HBL in large num- 
bers at high fluxes, due to their UV/X-ray peaked 
SEDs, but deeper X-ray samp les should reveal an increas- 
ingly large fraction of LBL (Gi ommi fc Padovanilll99"^ 
[Padova ni fc G iommi 1995a). In the second case the sit- 
uation is reversed, with the fraction of HBL expected to 
increase at lower fluxes i n the radio band an d be basically 
constant in the X-rays (jFossati et al.lll997[ ). 

Note that the question of which BL Lac class (HBL or 
LBL) is most numerous is not simply a demographical 
question but has also strong implications on the physics 
of relativistic jets. The frequency at which most of the 
synchrotron power is emitted, J/peak, in fact, shows a very 
large range in BL Lacs. Furthermore, Vpcak (x 7poak'^^' 
where 7peak is the Lorentz factor of the electrons emit- 
ting most of the radiation, 6 is the Doppler factor, and 
B is the magnetic field. Evidence towards a predomi- 
nance of LBL or HBL would then point towards Nature 
preferentially making jets which peak at IR/optical or 
UV/X-ray energies, thereby constraining also the physi- 
cal parameters of the jets. 

In order to test the two competing scenarios we have 
also derived the number counts for HBL. HBL should 
be defined in terms of the position of the frequency 
at which most of the synchrotron power is emitted, 
f^pcak- This requires multifrequency data and fitting 
their SEDs. One alternative is to use the X-ray-to- 
radio flux ratio, or the effective radio- X-ray spectral in- 
dex Qfrx- Bo na fide HBL, in fact, wer e shown to have 
Q!rx < 0.78 (P adovani fc Giommil 1 1996^ . iPadovani et al.l 
(|2003[ ). however, have found that a definition based solely 
on was not optimal, while the position of the sources 
on the aoxj ctro plane, which means using two (instead 
of one) effective spectral indices, appeared to be more 
sensitive to the synchrotron peak frequency, especially 
for FSRQ. They then defined a so-called "HBL box", 
derived b y using all HBL in the mult i -frequ ency AGN 
catalog of IPadovani. Giommi. fc Fiord (jl997l ). i.e., a re- 
gion of this plane within 2cr from the mean aro, ctox, and 
Qirx values of HBL. For comparison with previous results 
we adopt in this paper both definitions. 

The situation in the radio band is illustrated in Fig. [H 
The fraction of BL Lacs with a^x < 0.78 is 14lg%, while 
that in the HBL box is 9l|%. These values are slightly 
smaller than those given in IPadovani et al.l (|2003D due 
to our somewhat different BL Lac definition and also 
because in that paper we did not correct for the miss- 
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log (0.3 — 3.5 keV) [erg cm ^ s '] 

Fig. 5. — T he integral X-ray number counts for BL Lacs, adapted 
from (Padovan i & Giommi (1995 aj). Data for five X -ray selected 
samples are shown (see Padovani fc Giommil I1995al . for details) . 
Filled triangles represent the bivariate X-ray counts for the 1 Jy 
LBL with /x 3x 10"^'^ erg cm"'^ . The dotted line represents 
the DXRBS LBL with /x > 1.6 X 10"" erg cm-2 s"!, with filled 
squares showing the values at selected fluxes to show the errors 
involved. In both cases these define complete, X-ray flux limited 
LBL samples. The so lid line repre sents the X^-ray number counts 
for LBL predicted by [Giommi Padovani ( 1994) and revised by 
IPadovani fc Giomm I lll995al ). while th e dashe d line shows the pre- 
dictions of the blazar sequence llFbssati et al.lll997D . See text for 
details. 

ing corner in the radio flux - X-ray flux plane discussed 
above and in § 13.3! (Note that no correction is required 
for HBL sources, as their X-ray fluxes are such that they 
are all well away from that corner: see Fig. |3|). Interest- 
ingly, the fraction of HBL, according to both definitions, 
appears to be constant with radio flux, within the er- 
rors, at least where we have enough statistics (/r < 0.1 
Jy). Said differently, the ratio LBL/HBL is constant and 
~ 6 (using the ctrx definition, for consistency with pre- 
vious work). At t he radio fiuxes reached by DXRBS, 
iFossati et al.l ()1997( ) predicted a value ~ 2.1 (extrapolat- 
ing from their Tab. 3), that is a factor of 3 smaller. Fig. 
|6|(see § I4.3|) paints a strikingly similar picture for FSRQ. 

To see what happens in the X-ray band we have defined 
an X-ray flux limited sample of LBL by applying a cut 
to our DXRBS sample at /x ~ 1.6 x 10~^^ erg cm~^ 
s~^ (dot-dashed line in Fig. O bottom). The resulting 
11 objects form a complete. X-ray selected LBL sample 
which reaches approximately the same limiting flux as 
the EMSS sample. The X-ray number counts for this 
sample are shown in Fig. [5] (dotted line, filled squares) 
to be in extremely good agreem ent with the predictions 
of IPadovani fc Giommil (|1995af ) (solid line)^. Note also 
how the LBL/HBL ratio increases roughly sevenfold at 

^ These number counts p redictions were originally published by 
IGiommi fc Pad ovanil II 1994I) and then revis ed, but not published in 
this form. bv lPadovam fc Giommil Ill995al ). 
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Fig. 6.— The integral number counts for DXRBS FSRQ at 5 
GHz (6 cm). The soUd Une represents the total number counts, 
while filled circles show the values at selected fluxes to show the 
errors involved. The dotted-dashed line shows the counts without 
the correction which we appli ed t o make DXRBS equivalent to 
a radio flux-limited sample (§ 13.31 1. The open square represents 
the surface density for the 1 Jy sample while the open triangle 
represents the surface density for the PKS 0.25 Jy sample. The 
dashed line and filled triangles are the number counts for the FSRQ 
with Orx < 0.78, while the dotted line and filled squares are the 
number counts for the FSRQ in the HBL box. See text for details. 
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fainter fluxes, going from ^ 1 — 2% at /x ^ 10 
10"" erg cm-2 g"! t o ^ 10% at fx ^ 1 0"^'^ erg cm-^ 
s^^. For comparison. iFossati et al](|1997l ) predicted the 
opposite behavior, that is a slight (5 — 10%) decrease 
over the same flux range (dashed line in Fig. We also 
evaluated the X-ray LF for the DRXBS LBL sample and, 
again, found it to be in very go od agreem ent wit h the 
predictions made by iPadovaiiifc Giomni] (|1995aD (see 
their Fig. 4). 

All available evidence seems then to favor a scenario 
where HBL represent a small (~ 10%), constant frac- 
tion of the BL Lac population, contrary to the predic- 
tions of the blazar sequence. This adds to the recent 
quest ioning of the general validity of the blazar sequence 
'e.g. ,lGiommi. Menna. &: Padoyanilll999l:lPadovani et al 



200a ICaccianiea fc Marchal l2004 lAnton fc Browne 
2005' 'Nieppol a et all 120061 : iLandt et all I2006L see also 
Padovani (200^)" Note that in this case, unlike that 
for FSRQ and BL Lacs 16. 3p . relative number counts 
do translate into relative space densities, as none of the 
two classes evolve. 

4.3. FSRQ 

Fig. [6]presents the integral number counts at 5 GHz for 
the DXRBS FSRQ (solid line), compared to the values 
derived from the 1 Jy (Stickcl ct al. 1994, open square) 
and the PKS 0.25 Jy (jWall et all I2005i . open triangle) 
samples. The latter value has been derived by convert- 
ing the numbers given in iWall et al.l ((2005h to 5 Ghz 



assuming a,- ^ and by multiplying them by the ratio 
of DXRBS FSRQ with < 0.5 and those wi th a, < 0.4 
^ 1.2 ), to compensate for the fact that the IWall et al.l 
2005f ) sample includes only sources with < 0.4. As 
before, the filled circles show the values at selected bins 
to show the errors involved without crowding the plot. 

The following points can be made regarding Fig. [S| 
1. our number counts agree with previous estimates at 
250 mjy; they appear to be higher than previous surveys 
at IJy but, given the error bars, not significantly so (al- 
though this might be related to the effect discussed in 
S 13. 6p : 2. the number ratio between FSRQ and BL Lacs 
appear to be independent of radio flux and ~ 8 (tak- 
ing into account the possible excess of DRXBS FSRQ at 
high fluxes). As for BL Lacs, the southern and northern 
counts agree within the errors down to ^ 150 — 200 mJy, 
where the effect of the higher flux limit of the northern 
hemisphere starts to become relevant. 

We have also derived number counts for "HBL- like" 
FSRQ, or HFSRQ, defined as for BL Lacs in terms 
of their arx and their position in the aox, aro plane. 
These are shown in Fig. [6l The fraction of FSRQ with 
arx < 0.78 is 12^2%, while that in the HBL box is is 
4t^%. The fraction of HFSRQ, according to both defi- 
nitions, appears also to be roughly constant with radio 
flux, within the errors, as was the case for HBL. 

We note that our derived number counts have im- 
portant implications also for the study of the Cosmic 
Microwave Background (CMB), as DXRBS reaches ra- 
dio fluxes which are much fainter than those of the 
foreground sources detected in the WMAP catalogue. 
iGiommi et al.l (|2006l ) have indeed shown that our surface 
densities imply that a significant number of faint blazars 
are expected to contaminate CMB fluctuation maps as 
foreground sources, with important implications for the 
CMB fluctuation spectrum. 

5. EVOLUTION 

We study the evolutionary properties of DXRBS 
through the 14 /14 test (|Avni fc BahcalU 119801 : 
.Morris et al.l Il99ll ). a variation of the F/l^iax test 
(jschmidtl ligelT" that is the ratio between enclosed 
and available volume. Values of (14/14) significantly 
different from 0.5 indicate evolution, which will be 
positive (i.e., sources were more luminous and/or more 
numerous in the past) for values > 0.5, and negative 
(i.e., sources were less luminous and/or less numerous 
in the past) for values < 0.5. Moreover, one can also 
fit an evolutionary model to the sample by finding the 
evolutionary parameter which makes (14/14) = 0.5. 

We have computed 14/K, values for our sources taking 
into account our flux limits (at 6 and 20 cm and in the X- 
ray band) and the appropriate sky coverage. Statistical 
errors are given bycr = l/\/l2N ()Avni fc BahcaTl[l980l ). 
To have a first, simple estimate of the sample evolution 
we have also derived the best fit parameter r assuming a 
pure luminosity evolution of the type normally used, i.e., 
P{z) = P(0)exp[r(z)/T], where T{z) is the look-back 
time (the smaller r the stronger the evolution). 

Table [2] gives the sub-sample in column (1), the num- 
ber of sources in column (2), the mean redshift in column 
(3), (14/14) and r in columns (4) and (5) for our A cos- 
mology, and (14/14) and r in columns (6) and (7) for the 
empty Universe cosmology, for comparison with previous 
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TABLE 2 
DXRBS Evolutionary Properties. 









Ho = 70, = 


0.3, Ca = 0.7 


Ho = 50, 


= 0, Ha = 


Sample 


N 




(Vc/Va> 


r 




T 


FSRQ 


129 


J_ . ^ u. uo 




,-, r,7+0.05 

^•^'-0.03 






t SKQ, /6cm > lUU mjy 


102 


i.40 ± U.O/ 


O.bOo ± O.U2y 


n oQ+0.09 

'-'■■'•^-0.06 


U.Oyo ± 0.029 


n qr, + 0.13 

"•'-'^-0.06 


FSRQ, /x//r > 4 X 10-" 


101 


1.55 ±0.07 


0.724 ±0.029 




0.723 ± 0.029 


n -,c+0.02 

0-16_o.oi 


FSRQ, z < 1 


49 


0.71 ± 0.03 


0.629 ± 0.041 




0.640 ± 0.041 


n -14+0-05 

'^•^^-0.03 


FSRQ, Orx < 0.78 (HFSRQ) 


37 


1.79 ±0.16 


0.735 ± 0.047 


^■^^-0.03 


0.726 ± 0.047 




FSRQ, in HBL box 


15 


1.59 ±0.23 


0.740 ±0.075 


1Q+0.05 


0.724 ±0.075 


16+°°5 


BL Lacs 


24 


0.26 ±0.04=^ 


0.54 ±0.06'' 


c 


0.57 ±0.061" 


c 


BL Lacs, /scm > 100 mJy 


17 


0.29 ±0.05=^ 


0.42 ± 0.07'' 


c 


0.45 ± 0.071" 


c 


BL Lacs, /x//r > 4 X 10"" 


22 


0.26 ± 0.04=^ 


0.54 ±0.06'' 


c 


0.56 ±0.061" 


c 


BL Lacs, «rx < 0.78 (HBL) 


7 


0.35 ±0.05" 


0.54 ±0.10'' 


c 


0.56 ±0.10'' 


c 


BL Lacs, arx > 0.78 (LBL) 


17 


0.19 ±0.03^ 


0.54 ±0.07'' 


c 


0.57 ±0.07'' 


c 


BL Lacs, in HBL box 


6 


0.37 ± 0. 06^5 


0.48 ±0.13'" 


c 


0.50 ±0.13'' 


c 



Excluding the 7 sources without redshift Assuming z ■ 
"om 0.5: no evolution assumed*' Excluding 6 sources 
without redshift^ Excluding 1 source without redshift 



(z) for the 



sources 



without redshift^ (V'c/Va) not significantly different 



from 0.5: no evolution assumed*' Excluding 6 sources without redshift*' Excluding 2 sources without redshiffl^ Excluding 5 sources 



results. Note that (z) is calculated taking into account 
the effect of the sky coverage (see discussion in Paper II) . 
The main results are the following: 

1. DXRBS FSRQ evolve at the 5.4ct (4.8cr for an 
empty Ui ii verse cosmology) level, a well known re- 
sult ( e.g., [Padovani fc Urrvl[l99^IUrrv fc Padovanil 
119951 ). Their evolutionary parameter for the sim- 
ple case of pure luminosit y evolution is fully consis - 
tent with that derived by l Urrv fc Padovanil ()1995f ) 
for the 2 Jy FSRQ sample for an empty Universe 
(r = 0.23;HD; 

2. DXRBS BL Lacs do not evolve, i.e., their (V^/Va) 
value is not significantly different from 0.5 (and 
consequently r > 1). The results for BL Lacs 
are more uncertain because of the smaller num- 
ber statistics. The fact that ~ 30% of them have 
no redshift is less of a problem, as redshift affects 
Vc/Vi^ values way less than flux (a value equal to 
(z) was assumed in this case); 

3. The typical redshifts for FSRQ and BL Lacs are 
markedly different, with the former covering the 
0.2-4.7 range and having (z) = 1.82±0.08 and the 
latter having 0.04 <z< 0.73 and (z) = 0.26±0.04; 

4. The (K/Va) values for HBL and LBL (using both 
definitions; see § |4]) are not significantly different. 
This is a new result, which contradicts the com- 
monly accepted point of view that HBL and LBL 
have different evolutionary properties. Notice that 
for the first time we can study the evolution of 
HBL and LBL within the same sample. Previ- 
ous comparisons had been typically made between 
the 1 Jy (radio-selected) an d the EMSS sam ples 
(X-ray-selected) , although R ector et al.l ()2000f ) did 
study the dependency of (T4/14) on X-ray-to-radio 
flux ratios for the EMSS sample. Admittedly, the 
errors on the (V^/Va) values are rather large but 
this is the best that can be done at present; 

5. HFSRQ have mean redshifts and evolutions similar 
to those of the main FSRQ sample. Although the 



{Ve/Va) values appear shghtly larger for HFSRQ 
(but still within < 2a), this difference can be easily 
explained by the fact that the HFSRQ sample is 
free from the effect discussed in § 13. 3[ unlike the 
full sample (see below), and completely identified, 
as all still to be observed sources have Ofrx > 0.78 
and are outside the HBL box. 

To check for the effect of the corner of the radio flux 
- X-ray flux plane "missed" by DXRBS (§ ES]), which 
is necessary only if one wants to compare our results to 
those of purely radio flux-limited samples, one cannot 
take the approach we used in § 2] The Vc/Va test, in 
fact, involves redshift, and so a simple correction to the 
observed number of sources with a given radio flux will 
not be sufhcient. We have taken two complementary ap- 
proaches to address this: 1. we evaluated (T4/T4) for 
both FSRQ and BL Lac samples applying a cut in radio 
flux at 100 mJy. By doing this the corner that DXRBS 
is missing in the radio-X-ray flux plane in Fig. [3] shrinks 
considerably, although we do lose a factor of two in ra- 
dio flux depth. Tab. [2] shows that (14/14) for these 
higher radio flux samples is consistent with that for the 
full samples within ~ 1.3cr; 2. we defined a sub-sample 
with complete coverage of the radio-X-ray flux plane for 
/x//r > 4 X 10"^^ erg cm"^ Jy~^ (to the right of 
the dotted line in Fig. [31 bottom; this translates into 
an additional constraint on the X-ray flux limit, which is 
bound to be > 4 X 10^^^ x (/r/Jy) erg cm^^ s^^). By 
doing this DXRBS has complete coverage of a somewhat 
restricted region of the radio-X-ray flux plane down to 
our nominal radio flux limit. Tab. [2] shows that (14/14) 
for these sub-samples is slightly larger than for the full 
samples but still within ^ 2a for our adopted cosmol- 
ogy. Both these results show that the double X-ray/radio 
selection discussed in § 13.31 has only a small effect on 
the derivation of the evolutionary properties of DXRBS. 
This makes sense as our X-ray flux limit, which is above 
the one which would guarantee complete coverage of the 
radio-X-ray flux plane, is properly taken into account 
when calculating the volumes (but this was not the case 
when deriving number counts). 
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Fig. 7.— The banded {V^/V^) statistic, ((Vi - Vo)/{V^ - Vo)) 
versus Zo for DXRBS FSRQ. The horizontal dashed Une indicates 
the value of 0.5 expected under the null hypothesis of no evolution. 
The (1 a) statistical error bars, given by c = l/vl2JV, where A'^ 
is the number of objects in the subsample with z > Zo, are shown 
for selected redshifts. 

We have also assessed the effect of the still unidentified 
sources as follows. We have assumed that all remaining 
10 sources with < 0.5 are FSRQ and all remaining 
6 sources with > 0.5 are BL Lacs. (Note that based 
on the relative fraction of the different DXRBS classes 
the expected numbers would be 9 and 2 for FSRQ and 
BL Lacs respectively: see § 12. 3p . We then assigned a 
redshift equal to (z) and then added these sources to the 
FSRQ and BL Lac samples. The resulting (Vc/Va) values 
increase only by < 0.6a as compared to the values given 
in Tab. [51 which shows that our results are quite stable 
against the addition of the still unidentified objects. 

Given the available FSRQ statistics we can move be- 
yond the assumption of pure luminosity evolution and 
study in detail possible redshift dependencies. In the 
case of a pure (luminosity or density) evolution model, 
that is under the assumption that the rate of change is 
independent of cosmic epoch, the best fit evolutionary 
paramete r has to be t he same at all redshifts (see, e.g., 
[della Coc a et al.lll992f ). It then follows that (Vc/Va) has 
also to be constant. The simplest test we can do is then 
to see if, and how, (14/14) changes with redshift. We 
then split our FSRQ sample in six redshift bins so that 
each bin contained roughly the same number of sources 
(22 for the first three bins and 21 for the last three). 
{Vc/Va} is roughly constant and ~ 0.6 between z ~ 0.3 
and 2, while above this redshift it drops to 0.47 (from 
0.65 at z ~ 0.3), a change which is significant at the 
2tT level. The sign of the evolution also changes, from 
strongly (2.5(t) positive at z 0.3 to consistent with no 
evolution at z > 2. 

An alternative, more common way to check for a cos- 
mic epoch-dependent evolution is through the so-called 



banded (T4/K) statistic, i.e., ((K - 14)/(K - Vo)), 
where Vo is the cos mological volume enclose d by a red- 
shift Zo (see, e.g.. iDunlop fc Peacockl [TgQOO . This al- 
lows the detection of any high-redshift, possibly nega- 
tive evolution by separating it by the well-known strong, 
positive, low-redshift evolution. This is shown in Fig. 
[71 The change is very strong, with a highly significant 
drop in (T4/14) with redshift, and evolution vanishing by 
z > 1.3. {Ve/Va}, in fact, goes from being > 0.5 at the 
5.4cr level at z '--^ to values which are < 0.5, although 
not significantly so, for z > 2. At higher redshifts, in 
fact, the lack of sources makes it very difficult to accrue 
more meaningful statistics from this test. This is an in- 
herent limitation of DXRBS, due to the fact that the 
faintest X-ray sources are only detected in a relatively 
small area of the sky (see Fig. [T]). Since, on average, 
these are also the sources with the highest redshifts, this 
explains our relatively large error bars at z > 2. These 
results, however, are consist e nt, for example, with thos e 
of Dunlop & Peacock (19901 iJarvis fc Rawlingj (|2000D . 
and Arshakian ct al. (2006) . Stronger evidence for a red- 
shift cut-off comes from the evolution of the FSRQ lumi- 
nosity function f? I6.2p . 

As regards BL Lacs, due to their smaller number, we 
did two simpler checks: first, we divided the sample 
with redshift information into two bins containing an 
equal number of sources, above and below z = 0.264; 
the {Ve/Va} values for the two samples are both consis- 
tent with 0.5 within 2a; second, we applied the banded 
(Vc/Va) statistics for Zq = (whole sample) and Zq = 
0.264; again, the two values are consistent with 0.5, this 
time within 1.3a. In the following we then assume no 
evolution at all redshifts for the BL Lac sample. 

6. LUMINOSITY FUNCTIONS 

As for the case of the evolutionary properties, the dou- 
ble X-ray/radio selection discussed in § 13.31 has only a 
very small effect on the derivation of the DXRBS LFs. 
In particular, by adopting the two approaches discussed 
in § [SI (i.e., applying a radio cut at 100 mJy and defining a 
sub-sample with fx/ fi > 4 x 10""'^^ erg cm~^ s""'^) there 
is little, if any, change in the BL Lac and FSRQ LFs, 
typically well within la at a given radio power. This 
is because the "missing" low radio flux sources would be 
spread over a range of powers, making the overall impact 
smaller. Moreover, the volumes are properly calculated 
using our X-ray flux limit (see § [5]). All our LFs are 
derived using the 1/Vmax (in our case 1/T4) technique 
(lSchmidtlll968l) . 

6.1. BL Lacs 

The LF of DXRBS BL Lacs is shown in Fig. E As dis- 
cussed above, based on the (T4/T4) analysis, no evolution 
is assumed. Therefore, the shown LF is supposed to be 
epoch-independent. We have assumed z = (z) ^ 0.26 for 
the 7/24 BL Lacs without redshift. 

The LF is well fitted by a single power law of the 
form (/i(i-'r) oc PjT^''. Varying the binning, the differ- 
ential slope is in the range 2.05 < Br < 2.23. For a bin 
size AlogP = 0.6, which is representative, a weighted 
least-squares fit yields (j){P,) cx p-2.i2±o.i6 (^2 _ ^ 2 for 
3 degrees of freedom) . The total number density of BL 
Lacs in the range 7 x lO^'^ — 6 x 10^^ W/Hz, derived in- 
dependently of bin size from the integral LF, is 840 ± 100 
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Fig. 8. — The radio luminosity function of DXRBS BL Lacs (fille d 
points). Error bars correspond to Icr Poisson errors IIGehrelslll986l) . 
The solid line is a weighted least-squares fit to the data. 

Gpc-3. 

We have checked how our assumption on the missing 
redshifts affects the LF determination in two ways: 1. 
we have excluded the 7 objects without redshift from our 
computation and multiplied the LF by 24/17. The re- 
sulting LF is consistent with the previous one within the 
errors, which shows that the overall shape of the LF is 
not strongly dependent on us assuming z = (z) for the 7 
BL Lacs without redshift; 2. the fact that these sources 
show a featureless continuum might suggest that their 
redshifts could be the highest in the sample. The lack 
of features, in fact, could be related to stronger beam- 
ing and, therefore, to (on average) higher luminosities. 
We have then obtained lower limits on the redshifts of 
these sources from their V magnitudes by making use of 
the fact that BL Lacs are hosted b y ellipticals of almost 
constant luminosity (e.g., Urry et al. 2000), which makes 
them adequate standard candles. By assuming a value 
for the ratio of the jet /galaxy flux when a BL Lac appears 
featureless, one can estimate the apparent magnitude of 
the host galaxy, which gives in turn a lower limit on 
the redshift of the BL Lac (as the jet/galaxy ratio could 
be higher). We have taken a conservative limit on the 
jet/gala xy ratio of a f e ature less BL Lac of one, based on 
Fig. 1 of lLandt et al.l ()2002[ ). and have used the relation 
between apparent V magnitude and redshift for lum inous 
ellipticals published by Browne & Marchl (Hill). The 
resulting lower limits span the range 0.11 — 0.64, with 
(z) ^ 0.28, very close to the value we ass umed. With 
a jet / galaxy ratio of 10, which based on iLandt et al.l 
(|2002D could be more appropriate for a featureless BL 
Lac, we obtain lower limits in the range 0.24 — 1.47, 
with (z) ~ 0.65. However, even assuming this value for 
the missing redshifts changes very little the resulting LF, 
which is consistent with the previous one well within the 
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Fig. 9. — The radio luminosity function of DXRBS BL Lacs 
(filled points) compared to the predictions of a beaming model 
based on the 1 Jy luminosity function and evolution (solid line, 
lUrrv fc Padovani 1995 ). The open s quares represent the 1 Jy lu- 
minosity functi on (Stickel et al.ll991h . Error bars correspond to Icr 
Poisson errors IIGehrelslll986ll . The beaming predictions and the 1 
Jy points have been corrected to take into account o ur somewhat 
different definition of a BL Lac. For consistency with lStickel et all 
(1991) the DXRBS luminosity function has been de-evolved to zero 
redshift and an Hq = 50 km s~^ Mpc~^, Cm = Oi a^nd = 
cosmology has been adopted. See text for details. 

errors. 

Fig. m shows the LF of DXRBS BL Lacs for an 
empty Universe cosmology (filled points), to compare 
it with previous determinations and the predictions of 
unified scheme s. The fi gure shows also the 1 Jy LF 
(open squares. iStickel et al. 1991), and the predictions 
of a beam ing model based on the 1 Jy LF and evolution 
(solid line. lUrrv fc Padovanill995[ ). These show what one 
should expect to find when reaching powers lower than 
those use d to constrain the LF at the high end. Note 
that since [Stickel et ahl (1991) had de-evolved their LF 
to zero redshift using their best fit value (r = 0.32), we 
have done the same for DXRBS BL Lacs, for which we 

get T = 0A5tlfs- 

A few interesting points can be made: 1. the DXRBS 
and 1 Jy LFs arc in very good agreement in the region 
of overlap, despite the factor ~ 20 difference in limiting 
flux; 2. the DXRBS LF reaches powers about one order 
of magnitude smaller that those reached by the 1 Jy LF, 
as expected given point n. 1; 3. the DXRBS LF is in good 
agreement with the predictions of unified schemes, which 
means that the unification of BL Lacs and Fanaroff- 
Riley (FR) type I radio galaxies seems to work also at 
low powers. Note th at while the 1 Jy LF was fitted by 
0(Pr) oc p-2-53±o.i5 (IStickel et al.lll99lD . the DXRBS LF 
is somewhat flatter, with (j){Pr) oc p^-2-3i±o.i8 ^jr^j, 
size AlogP — 0.6; xi ^ 2.7 for 4 degrees of freedom; as 
before, the slope does not change, well within the errors, 
for different bin sizes) . This is explained by the fact that 
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Fig. 10. — The differential radio luminosity function of DXRBS 
FSRQ in a P X (f>{P) form in six different redshift bins (includ- 
ing roughly the same number of sources): 0.2 — 0.65, 0.65 — 0.90, 
0.90 - 1.22, 1.22 - 1.47, 1.47 - 2.0, and 2.0 - 4.72. The dashed 
line indicates, for reference, the luminosity function in the lowest 
redshift bin, while the dotted line shows what the highest redshift 
LP should be in the case of pure luminosity evolution with t equal 
to the best fit parameter given in Tab. |2]for the wh ole FSRQ sam- 
ple. Error bars correspond to la Poisson errors jGehrels 1986). 
See text for details. 

the BL Lac LF predicted by unified schemes flattens out 
for Pr ^ 10^^ W/Hz, a region which is better sampled 
by DXRBS. For this cosmology the total number density 
of BL Lacs in the range 10^** — 6 x lO'^^ W/Hz, derived 
independently of bin size from the integral LF, is 310±35 
Gpc~^, to be compared with the value of 40 Gpc~'^ in 
the range 6 x lO^"* - 3 x 10^^ W/Hz for the 1 Jy LF. 

Two caveats are worth mentioning: 1. as done above, 
we have assumed z = (z) ~ 0.26 for the 7/24 BL Lacs 
without redshift but, as before, the overall shape of the 
LF is not strongly dependent on this assumption; 2. as 
mentioned above (§ |4]) the definition of a BL Lac for the 
1 Jy and DXRBS samples is somewhat different. Since 
2/3 (16/24) of the DXRBS BL Lacs fulfill the 1 Jy defini- 
tion, the 1 Jy points and the beaming predictions based 
on them were multiplied by the factor 3/2 to take this 
difference into account. 

6.2. FSRQ 

The case for FSRQ is more complex, as we know from 
the ((Vo - Vo)/(Va - V;)) analysis (§ EJ that the evolu- 
tionary parameter is epoch-dependent, and therefore we 
cannot simply de-evolve the global LF to zero redshift. 
We have then studied the LF evolution as a function of 
redshift. To be able to do so in a meaningful way retain- 
ing also a significant number of sources per redshift bin 
we have divided our sample in six bins so that each bin 
contains roughly the same number of sources (22 for the 
first three bins and 21 for the last three). Finally, we 
have computed both the differential and integral LFs, as 



log (5 GHz) [W Hz"'] 

Pig. 11. — The integral radio luminosity function of DXRBS 
PSRQ in six different redshift bins (including roughly the same 
number of sources): 0.2 - 0.65, 0.65 - 0.90, 0.90 - 1.22, 1.22 - 1.47, 
1.47— 2.0, and 2.0 — 4.72. The dashed line indicates, for reference, 
the luminosity function in the lowest redshift bin, while the dotted 
line shows what the highest redshift LF should be in the case of 
pure luminosity evolution with r equal to the best fit parameter 
given in Tab. [2] for the whole FSRQ sample. See text for details. 



they give complementary information. 

Fig. [TO] shows the differential LF for DXRBS FSRQ 
in a P X (t){P) form. This is equivalent to the (f){MB) 
form normally used in the optical band and allows an 
easy separation of luminosity and density evolution as 
the former would simply translate the LF to the right 
(higher powers) with no change in the ordinate (num- 
ber), while the opposite would be true for the latter. 
The figure also shows (dotted line) what the highest red- 
shift LF should be in the case of pure luminosity evo- 
lution with T equal to the best fit parameter given in 
Tab. [2] for the whole FSRQ sample. Keeping in mind 
that for a single power law LF one cannot distinguish be- 
tween luminosity and density evolution, we assume that 
some luminosity evoluti on takes place, bas ed on studies 
in other bands (see, e.g.. lCroom et al.ll2004D . Despite our 
somewhat limited statistics, and related relatively large 
scatter, which prevent us from a more quantitative anal- 
ysis, three things are apparent: 1. most of the luminosity 
evolution happens at relatively low redshift, as already 
by z ~ 1.3 the increase in power is almost as high as 
expected in the highest redshift bin for a simple pure 
luminosity evolution model. Beyond z w 1 not much ac- 
tion seems to take place; 2. all but one of the six LFs can 
be fitted by a single power-law with slope in the range 
1.6 - 1.9 (0(Pr) oc P~^-; 4/6 are actually in the 1.6 - 
1.75 range), the exception being the highest redshift one, 
whose anomalous steepness is due to a single object. No 
trend of slope with redshift is present, within the errors; 
3. the total number of sources (maximum space density) 
appears to be roughly constant at low redshifts (z <1). 



16 



Padovani et al. 



o 

o 



A 
O 



1 1 1 . 


1 1 1 1 

/ % 




1 1 1 1 1 1 1 1 1 

- 

■ 

=^ ' ^ 














10=" 








lo"" 








• 1 






10='" 








1 


—3^ ' 





1 2 3 4 5 

redshift 

Fig. 12. — Integral number densities for DXRBS FSRQ as a 
function of redshift for four power ranges, from top to bottom: 

Pr > lO^S-S W/Hz, Pr > 10^7 W/Hz, Pr > lO^'^-S W/Hz, P,. > 

10^8 W/Hz. The dotted lines show the integral number densities 
expected in the case of pure luminosity evolution with t equal to 
the best fit parameter given in Tab. [2]for the whole FSRQ sample. 



Fig. [TT] shows the integral LFs. Here the sharp shift 
in power from z 0.4 to z ~ 1.3 is clearer, as is the sud- 
den stop at higher redshifts. There is also evidence, for 
10^^ W/Hz, of an increase in the number density 
of FSRQ with redshift (at a given power) and a hint of 
a decrease (with redshift) at larger powers. The evolu- 
tion of the number density is much better seen in Fig. [T^ 
which plots it as a function of redshift above four powers, 
labeled in the figure. The curves show an initial, strong 
increase at low redshift /powers, peaking at z ~ 1.5, fol- 
lowed by a decline at higher redshifts. Moreover, while 
at low powers the integral number densities are consis- 
tent with those expected in the case of pure luminosity 
evolution (dotted lines), at high powers (J' ^ 3 x 10^^ 
W/Hz) they are lower, pointing to a deficit of sources at 
high redshifts and luminosities as compared to the pure 
luminosity evolution scenario. All of the above suggests 
a high-redshift {z « 1.5) decline in the comoving space 
density of high-power (P > 10^^ W/Hz) FSRQ. These re- 
sults are similar to those obtained bv lWall et al.l (j2005l ) 
for the Parkes 0.25 Jy sample. 

We have also derived the LF of DXRBS FSRQ for an 
empty Universe cosmology, to compare it with previous 
determinations of the local LF of FSRQ and the pre- 
dictions of unified schemes. Given that the evolutionary 
parameter is epoch dependent (§[5]) and that therefore we 
cannot simply de-evolve the global LF to zero redshift, 
we have restricted ourselves to sources having z < 1. 
This is the redshift where the evolution appears to slow 
down considerably (see Figs. [TOland fTT|) and it gives us 
also a large enough number of objects. This sub-sample 
includes 49 sources and is characterized by r = 0.14;to 03 




log (5 GHz) [W Hz"^] 

Fig. 13.— The radio luminosity function of DXRBS FSRQ with 
z < I (filled points) de-evolved at zero redshift compared to the 
predictions of a beaming model based on the 2 Jy luminosity func- 
tion and evolution (solid line, Urry & Padovani 1995). The open 
squares represent the 2 Jy luminosity function. Error b ars rep- 
resen t the sum in quadrature of the la Poisson errors IGehrelsl 
119861) and the variations of the number density associated with a 
IfT change in the evolutionary parameter r. For consistency with 
lUrrv fc Padovanil ljl995,) an Hq = 50 km s~i Mpc"!, S7m = 0, and 
Qa = cosmology has been adopted. Both beaming model and 2 
Jy LF have been converted from 2.7 GHz assuming = 0. See 
text for details. 

(see Tab. Note that this evolution is as strong as 
found in the optical band for the 2dF QSO Redshift Sur- 
vey ( 2QZ) (r = 1/ki = 1/6.15 - 0.16: ICroom et al.l 
[2OOI . 

The resulting LF is presented in Fig. [T2] (filled points), 
which shows also the 2 Jy LF and the predictions of 
a beaming model based on the 2 Jy LF and evolution 
(solid line, lUrrv fc Padovanil 1 1995f ). These show what 
one should expect to find when reaching powers lower 
than those used to constrain the LF at the high end. 
The model and 2 Jy LF have been converted from 2.7 
GHz assuming ar = 0. 

A few interesting points can be made: 1. the 2 Jy and 
DXRBS LFs are in very good agreement in the region 
of overlap, despite the factor ~ 40 difference in limiting 
flux; 2. the DXRBS LF reaches powers more than one 
order of magnitude smaller that those reached by the 2 
Jy LF, as expected given point n. 1; 3. the DXRBS 
LF is in good agreeme nt with the p redictions of unified 
schemes of Urry & Pad ovanil (|1995f ) , although slightly 
above them in the first two luminosity bins, which means 
that the unification of FSRQ and FR II radio galaxies 
seems to work also at low powers. The observed differ- 
ences could be simply due to the fact that the predictions 
are based on a FR II LF derived from a very high flux 
(2 Jy) sample; 4. we are getting close to the limits of 
the FSRQ "Universe". In fact, as FSRQ are thought to 
be the beamed counterparts of high-power radio galax- 
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Fig. 14. — Top: the difTerential radio luminosity function in a 
P X (f>{P) form for BL Lacs (black) in two redshift bins (0 — 0.26; 
dashed and 0.26 — 0.73; dotted lines) and FSRQ in three redshift 
bins: 0.2 - 0.65 (red), 0.65 - 0.90 (green), 0.90 - 1.22 (blue) and 
at 2 = (magenta). The latter is based on the ^ < 1 sub-sample 
and was de-evolved to zero redshift using the appropriate evolu- 
tionary parameter (se e § a nd Tab. Error bars correspond 
to la Poisson errors IIGehrelslll986i ). Bottom: the integral radio 
luminosity function of BL Lacs in two redshift bins (dashed and 
dotted lines) and FSRQ in three redshift bins and at ^ = 0. 

ies, the low-luminosity part of their LF should end at 
relatively high powers. Assuming that the minimum lu- 
minosity inferred from the fit to the 2 Jy LF is correct 
(solid line in the figure, bas ed on the 2 Jy LF of FR II 
radio galaxies; see lUrrv fc Pa dovani 1995), then DXRBS 
is approaching that value. 

Note that the model prediction of a flattening of the 
FSRQ LF for < 10^^ W/Hz, a region which is better 
sampled by DXRBS, fits quite well the fact that the 2 Jy 
LF is 4){P,) cx p-2-3±0-3^ ^hile the DXRBS LF is flatter, 
with (t){P,) cx p-i-63±o.i6^ 

6.3. BL Lacs and FSRQ 

DXRBS allows us to compare the LFs of BL Lacs and 
FSRQ within the same sample, which has obvious advan- 
tages. Given the redshift dependence of the FSRQ LF 
(Fig. fTU]) and the fact that BL Lacs reach only z ^ 1, we 
do this in Fig. [T3]for the BL Lac sample, split into two 
equally large sub-samples at z = 0.26^, and the z < 1 
FSRQ. In the latter case, we show the three lowest red- 
shift bins of Fig. [TO] and the LF of the z < 1 sources 
de-evolved to zero redshift using the appropriate evolu- 
tionary parameter (see Tab. [5]). 

A few interesting points can be made: 1. the different 
evolutionary properties of FSRQ and BL Lacs are visu- 
ally apparent. Namely, splitting the BL Lac LF into two 

''' We have assumed that the sources without redshift are equally 
split between the two sub-samples and re-normalized the LFs ac- 
cordingly. 



redshift bins is equivalent to a simple luminosity split, 
with the two LFs over lapping without discontinuity (see 
iBeckmann et al.l l2003l for a similar result in the X-ray 
band). The FSRQ LFs, on the other hand, clearly dis- 
play an "evolution", with the LFs in different redshift 
bins being shifted with respect to one another, as al- 
ready discussed in 16.21 2. BL Lacs are ~ 50 times more 
numerous than FSRQ, the former having a total num- 
ber density - 1,100 Gpc'^ for P, > 7 x 10^3 W/Hz, 
the latter having a total number density ~ 23 Gpc~^ for 
Fr > 4 X lO^'^ W/Hz. This is not simply due to the fact 
that BL Lacs reach lower powers, as the BL Lac number 
density above 4 x lO^"' W/Hz is ~ 240 Gpc~^, i.e., a fac- 
tor ^ 10 larger than that of FSRQ above the same radio 
power, due to the fact that the FSRQ LF is flatter than 
that of BL Lacs (0(Pr) oc Pf^ '' vs. (/)(Pr) oc P,"^'^)- 
The fact that FSRQ are « 8 times more abundant in 
our sample (§ H]) is due to the fact that FSRQ evolve, 
which means they are more luminous and their fluxes are 
"boosted", and to our flux limit. BL Lacs are in fact ex- 
pected to "catch up" at lower radio flu xes, becoming the 
dominant blazar class below « 2 mJy (jPadovani fc Urrvl 
1992); 3. all of the above, including the rough overlap 
between FSRQ and BL Lac LF in the lO^^ < < lO^^ 
W/Hz regi me, is in perfect accordance with the unified 
schem es of lUrrv IFP adovani (1995) (see also Radovan^ 
[1991) . Namely, a scenario in which BL Lacs are beamed 
FR Is and FSRQ are beamed FR lis explains the larger 
number densities, lower radio powers, and steeper LF of 
the former (see also Figs. 191 and [TO)) . 

Finally, the observed LFs constrain a possible evo- 
lutionary link between the two classes. It h as in 
fact been suggested (e .g., iCavaliere fc D'Elial [20021 : 
[Bottcher fc Dermeii[200a) that a "genetic" link might be 
present between FSRQ and BL Lacs, with some of the 
former switching from a short-lived, high accretion rate 
regime to a long-lived, low accretion rate regime. This 
would imply a decrease in the FSRQ number density at 
lower redshifts, which is not seen and goes also against 
the evidence discussed in § [5] (see Fig. [7]) . Fig. [T3] shows 
in fact that the FSRQ LF evolves smoothly to z = 0, 
while that of BL Lacs has no redshift dependence. One 
could then infer that any evolutionary connection be- 
tween FSRQ and BL Lacs has then to be limited and 
cannot affect the bulk of the blazar population. One 
caveat, however, is that radio power makes up a very 
small fraction of the total, bolometric luminosity, and is 
also not likely to be related to the episodes of large accre - 
tion of cold gas discussed bv[Cavaliere fc D'Elia[ ([2002[ ) 
and [Bottcher fc Derme^ ([2003 1 

7. SUMMARY AND CONCLUSIONS 

We have used a well-defined, complete sample selected 
from the Deep X-ray Radio Blazar Survey (DXRBS) to 
probe the radio number counts, evolution, and luminos- 
ity functions of blazars down to 5 GHz fluxes (^ 50 mJy) 
and powers (^ lO^"* W/Hz) about one order of magnitude 
deeper than previously available. Our sample includes 
129 flat-spectrum radio quasars (FSRQ) and 24 BL Lacs 
detected in the radio and X-ray bands over < 2, 000 deg^. 
Great care has been taken in assessing the completeness 
of the sample and the effects of its double X-ray/radio 
selection, which are relevant for other on-going surveys 
as well. Our main results can be summarized as follows: 
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1. Our number counts agree with previous estimates 
at higher radio fluxes; the surface densities of BL 
Lacs and FSRQ reach ~ 0.06 deg~^ and ~ 0.6 
deg~^ respectively at /sghz ^ 50 mJy. 

2. The two blazar sub-classes have diflFerent evolution- 
ary properties. FSRQ evolve as strongly as opti- 
cally selected quasars, at least up to rcdshift w 1.5, 
with evidence of no evolution and a decline in space 
density at higher redshifts for high-power sources. 
BL Lacs, on the other hand, do not evolve. This is 
true also for high-energy peaked BL Lacs (HBL), 
at variance with some previous results based on X- 
ray selected samples, which had found evidence of 
negative evolution. FSRQ and BL Lac redshifts are 
also widely different, with the former covering the 
0.2-4.7 range with [z) = 1.82 ±0.08 and the latter 
having 0.04 <z< 0.73 and (z) = 0.26 ± 0.04. 

3. The observed radio luminosity functions are in 
good agreement with the predictions of unified 
schemes, with FSRQ getting close to their expected 
minimum power (w 5xlO'^y{Ho/50)^ W/Hz). The 
idea that blazars are radio galaxies seen with their 
jets close to our line of sight seems then to work 
also at relatively low powers. 

4. Despite the fact that the large majority of DXRBS 
blazars are FSRQ, BL Lacs are intrinsically ~ 50 
times more numerous, again in agreement with uni- 
fied schemes. 

5. The observed relative numbers of HBL and LBL 
are different from those predicted by the so-called 
"blazar sequence". The LBL/HBL ratio is roughly 



constant and 6 in the radio band, instead of the 
predicted decrease down to ^ 2 at our flux lim- 
its. The opposite behavior is seen in the X-ray 
band, where a marked (sevenfold) increase in the 
LBL/HBL ratio is observed going to lower fluxes, 
instead of the slight decrease expected. The avail- 
able evidence supports a scenario in which HBL 
are intrinsically a small minority (~ 10%) of all 
BL Lacs. 
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